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MEASURE-VALUED MASS EVOLUTION PROBLEMS WITH FLUX 
BOUNDARY CONDITIONS AND SOLUTION-DEPENDENT VELOCITIES 

JOEP H.M. EVERS*, SANDER C. HILLEt, AND ADRIAN MUNTEAN* 

Abstract. In this paper we prove well-posedness for a measure-valued continuity equation with solution- 
dependent velocity and flux boundary conditions, posed on a bounded one-dimensional domain. We generalize the 
results of [EHM15a] to settings where the dynamics are driven by interactions. In a forward-Euler-like approach, 
we construct a time-discretized version of the original problem and employ the results of [EHM15a] as a building 
block within each subinterval. A limit solution is obtained as the mesh size of the time discretization goes to zero. 
Moreover, the limit is independent of the specific way of partitioning the time interval [0, T]. 

This paper is partially based on results presented in [Evel5, Chapter 5], while a number of issues that were still 
open there, are now resolved. 

Key words. Measure-valued equations, nonlinearities, time discretization, flux boundary condition, mild 
solutions, particle systems 

AMS subject classifications. 28A33, 34A12, 45D05, 35F16 


1. Introduction. A considerable amount of recent mathematical literature has been devoted 
to evolution equations formulated in terms of measures. Such equations are used to describe sys¬ 
tems that occur in e.g. biology (animal aggregations [CFRT10, CCR11], crowds of pedestrians 
[CPT14], structured populations [DG05, GLMC10, CCGU12, AI05]) and material science (defects 
in metallic crystals [vMM14]). Many interesting and relevant scenarios take place in bounded do¬ 
mains. Apart from the examples mentioned above, these include intracellular transport processes, 
cf. [EHM15b, Section 1], and also manufacturing chains [GHS+14]. However, most works that 
deal with well-posedness of measure-valued equations and properties of their solutions treat these 
equations in the full space, see for instance also [BGCG06, CDF+11, TF11, CLM13, CCS15]. The 
present work explicitly focuses on bounded domains and the challenge of defining mathematically 
and physically ‘correct’ boundary conditions. 

In [EHM15a], we derived boundary conditions for a one-dimensional measure-valued transport 
equation on the unit interval [0,1] with prescribed velocity field v. A short-hand notation for this 
equation is: 


d , d 


( 1 . 1 ) 


We focused on the well-posedness of this equation, in the sense of mild solutions , and the con¬ 
vergence of solutions corresponding to a sequence (f n )ne N in the right-hand side. Some specific 
choices for (/ n ) n eN represent for instance effects in a boundary layer that approximate, asn-> oo, 
sink or source effects localized on the boundary. The boundary layer corresponds to the regions 
in [0,1] where the functions f n are nonzero. 


There are several reasons why we consider mild solutions rather than weak solutions. First of 
all, the mild formulation in terms of the variation of constants formula - see (2.19) - follows 
directly from a probabilistic interpretation, as was shown in [EHM15a, Section 6]. Therefore the 
choice for mild solutions is justified by a modelling argument. Secondly, usually uniqueness of 
weak solutions cannot be expected to hold, while mild solutions are unique when the perturbation 
(/i i—> / • /i) is Lipschitz. In [EHM15a], where the perturbation even has discontinuities, we still 
obtain uniqueness of the mild solution. This is one of the main results of [EHM15a]. In the works 
[AI05, GLMC10, CCGU12, GJMC12] a specific weak solution is constructed that is precisely the 
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mild solution that we obtain by different means. Finally, there is a technical advantage of using 
mild solutions. Most of our estimates are in terms of the dual bounded Lipscliitz norm || • ||g L , 
that will be introduced in §2.1. Because test functions do not appear explicitly, our calculations 
are often simpler than when weak solutions are considered. Moreover, our estimates are in fact 
uniform over test functions in a bounded set. 


In the present work, we propose and investigate a procedure to generalize the former results 
to include velocity fields that depend on the solution itself. Such generalization makes it possible 
to model in a bounded domain the dynamics governed by interactions between the ‘particles’; in 
particular we will be concerned with interaction terms of convolution type that are given by a 
weighted average over the whole population. 


The results in this paper hold for a source-sink right-hand side that is based on a function / 
that is an element of the space BL([0,1]) of bounded Lipschitz functions on [0,1], In [EHM15a], 
we worked with / : [0,1] —> R that is piecewise bounded Lipschitz, though. Hence, here we are 
able to describe absorption in a boundary layer, but not yet absorption on the boundary alone. 
In the discussion section of this paper, see §5.1, we comment on the possibilities to extend our 
results to / that is piecewise bounded Lipschitz. 


We consider (1.1) for velocity fields that are no longer fixed elements of BL([0,1]). Instead of 
v, we write v[p] for the velocity field that depends functionally on the measure p. The transport 
equation on [0,1] becomes 




(1.2) 


The aim of this paper is to ensure the well-posedness of (1.2), in a suitable sense. Because (1.2) is a 
nonlinear equation, establishing well-posedness is not straightforward. Here, we employ a forward- 
Euler-like approach that builds on the fundamentals constructed in [EHM15a]. We partition the 
time interval [0,T] and fix the velocity on each subinterval. That is, restricted to a subinterval, 
the velocity depends only on the spatial variable and not on the solution measure. Within each 
subinterval the measure-valued solution evolves according to the fixed velocity and the evolution 
fits in the framework set in [EHM15a]. A more detailed description of our approach is given in 
§3. We decrease the mesh size in the partition of [0,T] and estimate the difference between Euler 
approximations. The main result of this paper is the fact that this procedure converges. 


A forward-Euler scheme similar to ours is used in [PR13] for measures absolutely continuous 
with respect to the Lebesgue measure. Their results are extended to general measures in [CPT14, 
Chapter 7]. The difference between their work and ours is twofold: they use the Wasserstein 
distance and they work in unbounded domains. 

The references that directly inspired us are [CG09, Hool3, GLMC10]. The approach presented 
in this paper deviates from [Hool3], since we restrict ourselves to evolution on the interval [0,1], 
while [Hool3] considers [0,oo). Furthermore, our regularity conditions on the velocity - given in 
Assumption 3.1 - are weaker than in [Hool3]; cf. Remark 3.3. Moreover, [Hool3] restricts to veloc¬ 
ity fields that point inwards at 0. In this way, no mass is allowed to flow out of the domain [0, oo). 
In our approach, the fact that the flow is stopped at the boundary is encoded in the semigroup 
(Pt) t^O; irrespective of the sign of the velocity there; cf. §2.2. We consider it too restrictive to have 
a condition on the sign of the velocity at 0 or 1; in practice it is very difficult to make sure that 
such condition is satisfied when the velocity v[p] depends on the solution (like in e.g. Example 3.2). 


In this paper we limit our attention to a one-dimensional state space, [0,1], because in this case 
the (global) Lipschitz continuous dependence of the stopped flow on the time-invariant velocity 
field v is a rather straightforward property (see Section 2.2, Lemma 2.2). In higher dimensional 
(bounded) state spaces this seems much more delicate to establish. We comment on this in more 
detail after the proof of Lemma 2.2. One should note however, that the results on convergence of 
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the forward-Euler-like approach that we present do not depend on the dimensionality other than 
through the mentioned Lipschitzian property as presented in Lemma 2.2. 


This paper is organized as follows. Within each subinterval of the Euler approximation the 
dynamics are given by a fixed velocity, like in [EHM15a]. Therefore, we start in §2 by collecting 
the results of [EHM15a] that we require here: a number of properties of the semigroup (Pt)t^o 
and of the solution operator, called (Qt)t^o- The forward-Euler-like approach to construct solu¬ 
tions is introduced in §3, where we also state the main results of this paper: Theorems 3.10 and 
3.12, and Corollary 3.11. In plain words and combined into one pseudo-theorem, these results read: 


Theorem. The proposed forward-Euler-like approach converges as the mesh size of the time 
discretization goes to zero. The limit is independent of the specific way in which the time domain 
is partitioned. This approximation procedure yields existence and uniqueness of mild solutions to 
the nonlinear problem, and solutions depend continuously on initial data. 


A more precise formulation follows later. We prove these results in §4 using estimates between 
two Euler approximations of (1-1). In §5 we reflect on the achievements of this paper, discuss 
open issues and provide directions for further research. 

2. Preliminaries. This section contains a summary of the results obtained in [EHM15a] on 
which we shall build. Moreover, we mention the technical preliminaries needed for the arguments 
in this paper. 

2.1. Basics of measure theory. If 5 is a topological space, we denote by A4(S) the space 
of finite Borel measures on S and by A4 + (5') the convex cone of positive measures included in it. 
For x G S, S x denotes the Dirac measure at x. Let 

(mM : = J <t>dp (2.1) 

s 

denote the natural pairing between measures p G M(S) and bounded measurable functions (j). 
The push-forward or image measure of p under Borel measurable <f> : S —> S is the measure Qffp 
defined on Borel sets E C S by 

(<t>#p)(E):=p(<S>-\E)). (2.2) 

One easily verifies that (4>#//, <j>) = (p, <f> o $). 


We denote by C b (S) the Banach space of real-valued bounded continuous functions on S equipped 
with the supremum norm || • Hoc. The total variation norm || • ||tv on A4(S) is defined by 


|tv : = sup• 


|(/!,</>) (f G C b (S), |MI 





It follows immediately that for 4> : S —> S continuous, ||$#/x||tv ^ IIHItv- In our setting, S 
is a Polish space (separable, completely metrizable topological space; cf. [Dud04, p. 344]). It is 
well-established (cf. [Dud66, Dud74]) that in this case the weak topology on At (5) induced by 
C b {S) when restricted to the positive cone At + (5 I ) is metrizable by a metric derived from a norm, 
e.g. the Fortet-Mourier norm or the Dudley norm. The latter is also called the dual bounded 
Lipschitz norm , that we shall introduce now. To that end, let d be a metric on S that metrizes the 
topology, such that ( S , d) is separable and complete. Let BL(5', d) = BL(S) be the vector space 
of real-valued bounded Lipschitz functions on ( S,d ). For <j> G BL(S'), let 


Ml := sup 


d(x,y) 


x,y G S, x^y 


be its Lipschitz constant. Now 


IMIbl := Moo + Ml 


(2.3) 
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defines a norm on BL(5) for which this space is a Banach space [FM53, Dud66]. In fact, with this 
norm BL(S') is a Banach algebra for pointwise product of functions: 

110 • V'IIbl < ||0||bl II^IIbl- (2.4) 

Alternatively, one may define on BL(5') the equivalent norm 

II0I|fm := m ax(||0|| oo , |0| L ), 

where ‘FM’ stands for ‘Fortet-Mourier’ (see below). Let || ■ ||g L be the dual norm of || • ||bl on the 
dual space BL(5)*, i.e. for any x* £ BL(S')* its norm is given by 

II^Ibl := sup { | (x*,<f>) | |0£BL(S), ||0||bl<1}. 

The map y i-a with 7^(0) := (p,, 0) defines a linear embedding of At (S') into BL(S)*; see 
[Dud66, Lemma 6]. Thus II ’ IIBL induces a norm on A4 (S’), which is denoted by the same symbols. 
It is called the dual bounded Lipschitz norm or Dudley norm. Generally, ||/r||g L ^ II/^IItv for all 
y £ A4(5). For positive measures the two norms coincide: 

NIbl = M(S) = IHItv for all y £ A4~*~(S). (2.5) 

One may also consider the restriction to A4(5) of the dual norm II ’ IIfm of || • ||fm on BL(S)*. 
This yields an equivalent norm on A i(S) that is called the Fortet-Mourier norm (see e.g. [LMS02, 
ZahOO]): 


IImIIbl ^ IImIIfm ^ 2||/.i||g L . (2.6) 

This norm also satisfies ||/z||p M ^ IImIItv, so (2.5) holds for || • ||p M too. Moreover (cf. [HW09, 
Lemma 3.5]), for any x,y £ S, 

114 - 4IIbl = 2 + d^y) ^ = II 5 - - 4IIfm- (2-7) 

In general, the space Af(S) is not complete for || ■ ||g L . We denote by A4(S)bl its completion, 
viewed as closure of At (S) within BL(S)*. The space Af + (S) is complete for |j • ||g L , hence closed 
in M(S) and A4(S)bl- 

The || • ||g L -norm is convenient also for integration. In Appendix C of [EHM15a] some techni¬ 
cal results about integration of measure-valued maps were collected. These will also be used in 
this paper. The continuity of the map x S x : S —> A4 + (S)bl together with (C.2) in [EHM15a] 
yields the identity 


y = J 6 x dy(x) (2.8) 

s 

as Bochner integral in At (S)bl; for basic results on Bochner integration, the reader is referred to 
e.g. [DU77]. The observation (2.8) will essentially link ‘continuum’ (‘/it’) and particle description 
(‘ch,’) for our equation on [0,1]. 

2.2. Properties of the stopped flow. Let v £ BL([0, 1]) be fixed. We assume that a single 
particle (‘individual’) is moving in the domain [0,1] deterministically, described by the differential 
equation for its position x{t) at time t: 


x{t) = v(x(t)), 
z(0) = x 0 . 


(2.9) 
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A solution to (2.9) is unique, it exists for time up to reaching the boundary 0 or 1 and depends 
continuously on initial conditions. Let x( ■; Xo) be this solution and I xo be its maximal interval of 
existence. Define 


Td{Xo) ■= sup I xo G [0, oo], 

i.e. rg(xo) is the time at which the solution starting at Xo reaches the boundary (if it happens) 
when xo is an interior point. Note that rg(xo) = 0 when Xq is a boundary point where v points 
outwards, while rg(x o) > 0 when xq is a boundary point where v vanishes or points inwards. 


The individualistic stopped flow on [0,1] associated to v is the family of maps 4> t : [0,1] -A [0,1], 
t ^ 0, defined by 


$t(zo) 


x{t;x o), if t G I Xo , 

x(rg(xo); Xq), otherwise. 


( 2 . 10 ) 


To lift the dynamics to the space of measures, we define P t : A4([0,1]) —>- At([0,1]) by means 
of the push-forward under $ t : for all pi G A4([0,1]), 

P t p := $ t #p = p 0$-!; (2.11) 

see (2.2). Clearly, P t maps positive measures to positive measures and Pt is mass preserving on 
positive measures. Since the family of maps ( c E > t)t>o forms a semigroup, so do the maps Pt in 
the space A4([0,1]). That is, (Tt)t^o is a Markov semigroup on _Ad[0,1] (cf. [LMS02]). The basic 
estimate 


II-PismIItv < IImIItv 


( 2 . 12 ) 


holds for p G M ([0,1]). 

In the rest of this section we summarize those properties of (Pt)t^o that are needed in this paper. 
We first recall Lemma 2.2 from [EHM15a]: 

Lemma 2.1 (See [EHM15a, Lemma 2.2]). Let p G A4([0,1]) and t,s G R + . Then 
(i) \\P t p - PspWbl < IMI°o ||m||tv \t - s|. 

(™) II^mIIbl < max(l, |$ t | L ) \\p\\^ h ^ el ,, l Lt ||/z||| iL . 

To distinguish between the semigroups on A4([0,1]) associated to v,v' G BL([0,1]), we write 
P v and P v , respectively. Analogously, we distinguish between the semigroups (&t)t^o and 
($)’ ) t ^o on [0,1] and between the intervals of existence 1% and 1% associated to (2.9). 

Lemma 2.2. For all p G A4([0,1]), v, v' G BL([0,1]) and t G Mg' 

\\ p tT ~ -P^VIIbl < b ~ v'\\oot\\p\\TY e Lt , 


where L := min(|u|L , |t/|l)- 

Proof. For any if G BL([0,1]), we have 

| (</>, Pfp - pfp) I = I o - (j, o , p) K |</>|l ||$? - Iloo IImIItv, (2.13) 

hence 

\\p t v p-pfp\\^^ II^-^'IIocIHItv. 


Let x G [0,1]. 


(2.14) 
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Case 1: t e fl . 


!*?(*)-*? (*)l = 


^(x))-</($* (x))ds 


l 

< |v|l J\$ v s {x) - (a:) I ds + ||v - v'W^t. 


Gronwall’s Lemma yields 


l$?(z) | < \\v -v'Woote^*, 


(2.15) 


for all x £ [0,1]. Due to the symmetry of (2.15) in v and v', the same estimate (2.15) can be 
obtained with |?/|l instead of |x|l, and hence, we can write min(|'u|i i , |f / |l) in the exponent. This 
observation yields, together with (2.14), the statement of the lemma. 

Case 2: t qL If. We extend v : [0,1] -A R to v : R —> R by defining v(x) := v(0) if x < 0 and 
v(x) := x(l) if x > 1. Then v is a bounded Lipschitz extension of v such that ||n||oo = ||u||oo and 
|u| L = Ml- Let : R -A R be the solution semigroup associated to the unique (global) solution 
to (2.9) with v replaced by v and with initial condition to be taken from the whole of R. We 
extend v' analogously to v'. 

Irrespective of whether t £ 1% or t $. If , and whether in the latter case <&}(x) = 4>)’ (x) or 
$ £ (x) (x) , the following estimate holds 

!$?(*)-$?'(*)! <!$*(*)(2.16) 
for all x £ [0,1]. estimate (x) — (x)| using the same ideas as in (2.14) and (2.15) and obtain 

11-PiV - L^VIIbl < Ik — v'Woot ||/t||tv exp (min (I'D | l , MIl) t). 


The statement of the lemma follows from the equalities Ml = |v|l, |w / |l = MM |k — i/||oo = 
Ik — |joo and Equation (2.16). The case t $. 1% is analogous. □ 


Remark 2.3. The definition of stopped flow in state spaces of dimension two and higher and 
establishing elementary properties of its lift to measures is more delicate than the one-dimensional 
case presented above. Consider an open domain Ll C R”, n > 2 (with sufficiently smooth bound¬ 
ary). Let LI be its closure and let v £ BL(fi,R n ) be a velocity field on LI. Solutions to the initial 
value problem (2.9) with xq £ Ll still exist for some positive time, but in this higher dimensional 
setting it may happen that trajectories of the flow in Ll defined by v are partially contained in the 
boundary dLl or even only ‘touch’ dLl. So reaching the boundary in finite time is not equivalent to 
‘leaving the domain’. 

By means of the Metric Tietze Extension Theorem one can extend v to v £ BL(R",R”) with 
preservation of Lipschitz constant and supremum norm, e.g. component-wise. Let Xy{t\x o) be the 
corresponding solution starting at Xq at t = 0. The exit time or stopping time of the solution 
starting at Xq could then be defined as 

Tq{x o) := inf{f > 0 | x t -(t; xfi) £ R” \ Ll}, 

with the convention that the infimum of the empty set is set to +oo. One must show that this value 
is independent of the particular extension v that was chosen. rf(x o) now replaces the similarly 
denoted termination time of the solution that was used above for the one-dimensional case. The 
stopped flow can then be defined as in (2.10). 

Then for t > 0, 

t Ar| (x) 

$t(x) = x+ j v($ v s (x))ds, 

0 


(2.17) 
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where t/\Tg(x) denotes the minimum oft and Tg(x). The equivalent of Lemma 2.2 can be obtained 
from (2.17) by means of Gronwall’s Lemma essentially, once one knows that the stopping time 
satisfies for fixed x £ fl a Lipschitz estimate of the form 

1 1 A Tg(x) - t A Tg (z)| < C\\v - v'Woot. (2.18) 

We are not aware of results in the literature that provide estimates like (2.18). Neither did we 
succeed in establishing such an estimate ourselves. The rich possible dynamics of solutions of 
higher dimensional systems of non-linear differential equations may even make such global Lips¬ 
chitz dependence of the velocity field impossible in general. Thus, a generalization of this part of 
the paper to two and higher dimensional state spaces seems not straightforward. 

2.3. Properties of the solution for prescribed velocity. We consider mild solutions to 
(1.1), that are defined in the following sense: 

Definition 2.4 (See [EHM15a, Definition 2.4]). A measure-valued mild solution to the 
Cauchy-problem associated to (1.1) on [0,T] with initial value v £ A4([0,l]) is a continuous map 
H : [0, T] A4([0, 1])bl that is || • ||t y-bounded and that satisfies the variation of constants formula 

t 

Ii t = P t v + J Pt- s Ff(ns)ds for all t G [0,T]. (2.19) 

o 

Here, the perturbation map Ff : A4([0,1]) —> -M([0,1]) is given by Ff(p) := f ■ pt. 

We showed in [EHM15a] that mild solutions in the sense of Definition 2.4 exist, are unique and 
depend continuously on the initial data. We repeat those results in the following theorem. 

Theorem 2.5. Let f : [0,1] —> R be a piecewise bounded Lipschitz function such that v(x ) ^ 0 
at any point x of discontinuity of f. Then for each T ^ 0 and /io G A4([0,1]) there exists a unique 
continuous and locally || • ||t y-bounded solution to (2.19). Moreover, there exists Ct > 0 such that 
for all initial values G A4([0,1]) the corresponding mild solutions fi and p! satisfy 

IlMt - A^IIbl < C'rll/io - MoIIbl (2-20) 


for all t £ [0, T\. 

Proof. See [EHM15a, Propositions 3.1, 3.3 and 3.5] for details. □ 

In this paper, we restrict ourselves to those functions / that are bounded Lipschitz on [0,1]; 
see §5.1 for further discussion on the need of this restriction. Let v £ BL([0,1]) and / G BL([0,1]) 
be arbitrary. For all t ^ 0, we define Q t : A4([0,l]) —> A4([0,l]) to be the operator that maps 
the initial condition to the solution in the sense of Definition 2.4. Theorem 2.5 guarantees that 
this operator is well-defined and continuous for || • ||r T . Moreover, Q preserves positivity, due to 
[EHM15a, Corollary 3.4], 

In the rest of this section, we give an overview of the properties of the solution operator Q. 

Lemma 2.6 (Semigroup property). The set of operators (Qt)t^o satisfies the semigroup prop¬ 
erty. That is, 


Qt Qs iz — Qt+s [z 


for all s,t ^ 0 and for all p G A4([0,1]). 
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Proof. The proof follows the lines of argument of [S94, p. 283]. We consider 

t + S 

Qt-\-sft Qt Qs ft — Pt-\-sft ^ Pt-\-s—a Ff (Q(T ft) d(J 

0 

t 

- PtQsft- J Pt-a PfiQa Qs ft) da , (2.21) 

0 

and observe that 

S 

P t Qsft=PtP s ft + Pt J P S -a FfiQa ft) da 

0 

s 

= Tt+ s /i + J ’ Pt+s-a Ff(Q<j fj.) do. (2.22) 


Because / £ BL([0,1]), the map o P s -a Ff{Qa ff) is continuous and hence it is measurable. 
Therefore, the second equality in (2.22) holds due to [EHM15a, Equation (C.3)]. A combination 
of (2.21) and (2.22) yields that 


t-\-S t 

Qt-\-sft Qt Qs ft — ^ P t-\-s—c r Pf(Q<r ft) da J P t—cr Ff(Qcr Qs ft) da 

s 0 

t 

= JP t _ a (. F f {Q a+s n) - F f (Q a Q s /j,)) do. (2.23) 

o 

To obtain the last step in (2.23), we use the coordinate transformation r := o — s in the first 
integral and subsequently renamed the new variable r as a. We estimate the total variation norm 
of (2.23) in the following way: 

t 

||Qt+sM — Qt Qs m||tv ^ J II Pt-a (Ff(Qa+s n) — Ff(Q a Q s fj,)) ||tv do 

0 

t 

< K QsM) ||tv do 

o 

t 

^ Halloo J^ ||/?cr+s M QaQs A^||tV do. 

0 

Here, we used [EHM15a, Proposition C.2(iii)] (noting that the integrands are continuous with 
respect to o) in the first line, (2.12) in the second line and the fact that / £ BL([0,1]) C Cb([0,1]) 
in the last line. Gronwall’s Lemma now implies that \\Qt+ s ia — Qt Qs m||tv = 0 for all s, t ^ 0. □ 


Lemma 2.7. For all /i £ A1([0, 1]) and s,t ^ 0, we have that 

WQtH - Qsn\ Ibl < IImIItv -(11/1100 + Halloo) • e^l- max ^) • \t - s\. 


Proof. The statement of this lemma is part of the result of [EHM15a, Proposition 3.3] .□ 
Lemma 2.8. For all fj, £ A1([0, 1]) and t ^ 0, we have that 
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(i) IIQtPlkv < IIpIItv expdl/Hoo t), and 

(a) \\Qtp\\ B L < IIpIIbl exp(|v| L t + \\f\\ B Lte^ Lt ). 

Proof, (i): This estimate is given in [EHM15a, Proposition 3.3]. 

(ii): By applying [EHM15a, (C.l)] and Lemma 2.1 (ii) we obtain from (2.19) the estimate 

t 

II Qt pIIbl < ex p(|f |l t) ||p||bl + J exp(|n| L (t - s))||/||bl ||QsP||bl ds. 

o 

Gronwall’s Lemma now yields the statement of Part (ii) of the lemma. □ 

Corollary 2.9. For all p, ^ £ A4([0,1]) and t ^ 0, we have that 

WQt P Qt HIbl < IIp-HIbl ex p(|u| L t + II/Hbl t e^ L 4 ). 


Proof. Apply Part (ii) of Lemma 2.8 to the measure p — v £ A4([0,1]). □ 

We write Q v and Q v to distinguish between the semigroups Q on A4([0,1]) associated to v £ 
BL([0,1]) and v' £ BL([0,1]), respectively. 

Lemma 2.10. For all v, v' £ BL([0,1]), p £ M ([0,1]) andt^O, the following estimate holds: 

II QtP - Qt' H IIbl < \\v ~ f'Hoo IIpIItv exp(Lt + \\f\\ B Bte Lt ) ■ [t + t 2 |l/lloo e ll/l|o ° *]> 
where L := min(|?;|L , |i/|l)- 

Proof. We have 

t 

WQtV-Qt'vWBL < ||P t V-P/p|lBL + J\\ P t-s F f{Q V sP) - Pt- s F f {Q V s v)\\b B ds. (2.24) 

0 

Lemma 2.2 provides an appropriate estimate of the first term on the right-hand side. For the 
integrand in the second term, we have 

|| P?- s Ff(QM - P t v - s F f (Q v s ' h)\\bl < \\ p t-s P f(QsP) - Pt- s Ff{Q» IIbl 

+ || P?L a F f {QW - Pt s F f (Q V s p)||k 

^ ll v — u, ||oc (t — s) ||A 7 /(QsP)II tv e L> ' t ~ s ' > 

+ el-'l- IIF/(Q» - FfiQl'n )||* L , (2.25) 

due to Lemma 2.2 and Lemma 2.1(h). We proceed by estimating the right-hand side of (2.25) 
and obtain 

IIP^EHQM - P?s F f{Q V s p)\\bb < \\v - v'lloo (t a) ll/lloo IIpIItv ell'll- s e L 

+ e i«'iL(t-)|| / ||BL||g>-Q;'p||B L , (2-26) 

where we use Part (i) of Lemma 2.8 in the first term on the right-hand side. Since the estimate 
in (2.26) is symmetric in v and v' , we can replace |i/|l by L. 

Substitution of the result of Lemma 2.2 and (2.26) in (2.24) yields 

IIQtP - Q V tP IIbl < ll« - t/||oo * IIpIItv e Lt (1 +1 H/lloo e"' 11 ” *) 

t 

+ eLt ii/Hbl J iiq>-q:'pHb L da . 

0 

The statement of the lemma follows from Gronwall’s Lemma. □ 
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3. Measure-dependent velocity fields: main results. This section contains the main 
results of the present work. We generalize the assumptions on v from [EHM15a] in the following 
way to measure-dependent velocity fields: 

Assumption 3.1 (Assumptions on the measure-dependent velocity field). Assume that v : 
M([0,1]) X [0,1] ->K is a mapping such that: 

(i) v[p] £ BL([0,1]), for each p £ At([0,1]). 

Furthermore, assume that for any R > 0 there are constants Kr, Lr, Mr such that for all 
p, v £ AfQO,1]) satisfying ||/i||tv ^ R and, |M|tv ^ R, the following estimates hold: 

(ii) |M/u]||oo < k r , 

(Hi) | v[p] | L ^ L r , and 

M || v[p\ - uMIloo ^ Mr \\p - p||* L . 

Example 3.2. An example of a function v satisfying Assumption 3.1 is: 

v[p\(x) := [ JC(x - y) dp{y) = (1C*p)(x), (3.1) 

[ 0 , 1 ] 

for each p £ At([0,1]) and x £ [0,1], with 1C £ BL([—1,1]). This is a relevant choice, because it 
models interactions among individuals. 

Remark 3.3. Parts (ii) and (Hi) of Assumption 3.1 are an improvement compared to [HoolS]. 
There, the infinity norm and Lipschitz constant are assumed to hold uniformly for all p £ Af ([0,1]); 
cf. Assumption (FI) on [HoolS, p. )0[. We note that the convolution in Example 3.2 satisfies As¬ 
sumption 3.1, but does not satisfy Assumption (FI) in [Hool3[. They require a uniform Lipschitz 
constant because their Lemma f.3 is an estimate in the || ■ ||g L -norm for which Part (ii) of our 
Lemma 2.1 is used. Our counterpart of Lemma f.3 in [Hool3] is Lemma 3.). We give an estimate 
in terms of the || • ||t \-norm using (2.12) which does not involve the Lipschitz constant. 

Our aim is to prove well-posedness (in some sense yet to be defined) of (1.2). That is, 

on [0,1]. As said in §2.3, we restrict ourselves to / that is bounded Lipschitz on [0,1]. 

We now introduce the aforementioned forward-Euler-like approach to construct approximate so¬ 
lutions. Let T > 0 be given. Let N ^ 1 be fixed and define a set a C [0,T] as follows: 

a := {t.j £ [0, T] : 0 < j < IV, t 0 = 0, t N = T, tj < t j+1 }. (3.2) 

A set a of this form is called a partition of the interval [0, T] and N denotes the number of subin¬ 
tervals in a. 

Let p 0 €E M.{[ 0,1]) be fixed. For a given partition a := {t 0 ,... Rn} C [0,T], define a measure¬ 
valued trajectory p £ (7(10, T]; Af([0,1])) by 

{ Pt ■■= Qt-tj Th > (tj » tj+i ]; 

Vj ■■= v[p tj ]; (3.3) 

Pt -o = Mo, 

for all j £ {0,..., N — 1}. Here, (Qt)t^o denotes the semigroup introduced in §2.3 associated to 
an arbitrary v £ BL([0,1]). Note that by Assumption 3.1, Part (i), Vj = v[p tj ] £ BL([0,1]) for 
each j. 
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We call this a forward-Euler-like approach, because it is the analogon of the forward Euler method 
for ODEs (cf. e.g. [But03, Chapter 2]). Consider the ODE dx/dt = v{x ) on R for some (Lipschitz 
continuous) v : R -A R. The forward Euler method approximates the solution on some interval 
|_i] by evolving the approximate solution at time tj , named Xj , due to a constant velocity 
v{xj). That is, x(t) « Xj + (t — tj) ■ v{xj) for all t £ (tj,tj+ 1 ]. 

In (3.3), we introduce the approximation fit, where fit results from fit , by the evolution due to the 
constant velocity field v [fit, ] ■ The word constant here does not refer to v being the same for all 
x £ [0,1], but to the fact that v corresponding to the same fit, is used throughout (tj,tj+ 1 ]. 

The conditions in Parts (ii)-(iv) of Assumption 3.1 are only required to hold for measures in 
a TV-norm bounded set, in view of the following lemma: 

Lemma 3.4. Let fig £ _A4([0,1]) be given and let v : AI([0,1]) x [0,1] —> R satisfy Assumption 
3.1(i). For a given partition a := {to,... An} C [0, T], let fi £ C([0, T]; _M([0,1])) be defined by 
(3.3). Then the set of all timeslices of fi, that is 

A '■= {fit '■ t £ [0,T]}, 

is bounded in both || • ||tv and II • Hbl- The bounds are independent of the choice of a. 

Proof. Fix j £ {0,..., N — 1} and let t £ (tj,tj+ 1 ]. By Part (i) of Lemma 2.8, we have that 

11 Ah 11 TV = II Qt-tj MtjllTV < ll/h^TV ex P(ll/l|oo {t — tj)) 

< ll/h, 11 tv expdl/Hoo (t j+ 1 - tj)) 

for all t G (tj,tj+ 1 ]. Iteration of the right-hand side with respect to j yields 

j 

11 Ah 11 tv < IIa^oIItv Ifexpdl/Hoo^i+i - L)) = ||mo||tv expdl/Hoo {t j+1 - t 0 )). 

2—0 

Hence, for all t G [0,T] 

IIA 4 *IItv < IImoIItv expdl/Hoo (Ajv -t 0 )) = ||/z 0 ||tv expdl/Hoo T). 

This bound is in particular independent of t, N and the distribution of points within a. The 
bound in || • ||g L follows from the inequality |MIb L ^ IMItv that holds for all v £ A4([0,1]). □ 

In this paper we construct sequences of Euler approximations, each following from a sequence 
of partitions (ak)ke N that satisfies the following assumption: 

Assumption 3.5 (Assumptions on the sequence of partitions). Let (ak)ke n be a sequence of 
partitions of [0, T] and let (Nk)keN C N be the corresponding sequence such that each ak is of the 
form 

a k := {t k j £ [0, T] : 0 < j A N k , t k 0 = 0, t k Nk = T , t k < t k +1 }. (3.4) 

Define 

M := max t k +1 — t k (3.5) 

for all k £ N. Assume that the sequence {M^)k^n is nonincreasing and M^ —>0 as k —>• oo. 

Example 3.6. The following sequences of partitions satisfy Assumption 3.5: 

• For all k £ N, take Nk := 2 fc , and let t k := jT/2 k for all j £ {0,..., Nk}. This implies 
that = T/2 k for all k £ N. This specific sequence of partitions was used in [Evel5, 
Chapter 5]. 
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• Fix q £ N + . For all k £ N, take N\ := q k , and let t k := jT/q k for all j £ {0,..., Nk}. 
This implies that M^ = T/q k for all k £ N. In the discussion section of [Evel5, Chapter 
5], the results of the current paper were conjectured to hold for this case. 

• For all k £ N, take Nk := k+ 1, and let t k := jT/(k+ 1) for all j £ {0,. .., Nk}. This im¬ 
plies that M^ = T/(k + 1) for all k £ N. This is an elementary time discretization (with 
uniform mesh size) used frequently when proving the convergence of numerical methods. 

• Let «o be a possibly non-uniform partition of [0, T]. Construct the sequence (otk)ke n i n 
such a way that any ak+i is a refinement of ak- That is, ak +i C ak for all k £ N. 
Elements may be added in a non-uniform fashion to obtain ak +i from ak, as long as 
M—r 0 as k —> oo. In this case (Nk)ket* is automatically nondecreasing. 

Also, some less straightforward sequences of non-uniform partitions are admissible, in which subse¬ 
quent partitions are not refinements. See for example Figure 3.1, in which two subsequent elements 
from the sequence (afe)fcew are given. These elements could indeed occur, since M^ k+1 l < M^ k \ 
This example is rather counter-intuitive, as there is a local growth of the mesh size at the left-hand 
side of the interval [0,T] when we go from ak to ak+i- Note that even lV( fc+1 ) < . However, 

admissibility of a sequence of partitions is only determined by the local ordering of the maximum 
mesh spacing (i.e. the condition M^ k+1 ^ ^ M^) and its long-time behaviour: M^ —> 0 as 
k -A oo. 


Oik 



Olk+l 


o — 

7 




ifc +1 J .&+1 ifc +1 

z O l 2 


Figure 3.1. Two possible subsequent partitions in a sequence ( a.k)keN satisfying Assumption 3.5. 


Remark 3.7. Assumption 3.5 implies that Nk -A oo as k -a oo. 

If (M( fc ))k e N is not nonincreasing, but still -A 0, then it is possible to extract a subsequence 
such that is nonincreasing. 

We define a mild solution in this context as follows: 

Definition 3.8 (Mild solution of (1.2)). Let the space of continuous maps from [0,T] to 
jM([0, 1])bl be endowed with the metric defined for all p, v £ C([0, T]; Ad([0,1])) by 

sup \\pt-vt\\BL- (3-6) 

te[o,T] 

Let (cifc)feeN be a sequence of partitions satisfying Assumption 3.5. For each k £ N, let p k £ 
C([0, T]; Ad([0,1])) be defined by (3.3) with partition ak- Then, for any such sequence of parti¬ 
tions (ttfe)fcg n, any limit of a subsequence of (p k )keN is called a (measure-valued) mild solution of 
( 1 . 2 ). 

The name mild solutions is appropriate, because they are constructed from piecewise mild so¬ 
lutions in the sense of Definition 2.4. 

Remark 3.9. Consider the solution of (3.3) for any partition a C [0,T]. Mass that has 
accumulated on the boundary can move back into the interior of the domain whenever the velocity 
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changes direction from one time interval to the next. This is due to the definition of the maximal 
interval of existence I xo and the hitting time tq{x q) in §2.2. 

In the rest of this paper we focus on positive measure-valued solutions, because these are the 
only physically relevant solutions in many applications. The main result of this paper is the fol¬ 
lowing theorem. 

Theorem 3.10. Let po £ M + ([0, 1]) be given and let v : A4([0,1]) x [0,1] —> R satisfy As- 
sumption 3.1. Endow the space C([0,T];-Ad([0,1])) with the metric defined by (3.6). Then, there 
is a unique element of C([0, T]; A4 + ([0,1])) with initial condition po, that is a mild solution in the 
sense of Definition 3.8. That is, for each sequence of partitions satisfying Assumption 

3.5, the corresponding sequence (p k )ken defined by (3.3) is a sequence in C([0, T\; A4 + ([0,1])) and 
has a unique limit as k —>■ oo. 

Moreover, this limit is independent of the choice of (ak)ken- 

Corollary 3.11 (Global existence and uniqueness). For each po £ A4 + ([0,1]) and v : 
_A/f([0,1]) x [0,1] — > R satisfying Assumption 3.1, a unique mild solution exists for all time t ^ 0. 

Theorem 3.12 (Continuous dependence on initial data). For all T > 0 and R > 0 there is a 
constant such that for all po,vo £ Al + ([0,1]) satisfying ||/xo||tv ^ R and ||^o||tv ^ R, the 
corresponding mild solutions p, v £ C([0, T]; A4 + ([0,1])) satisfy 

sup || p T — za-IIbl ^ Cft T Who ~ ^oIIbl- 

tG[0,T] 


The proofs of these theorems and this corollary are given in the next section, §4. The key idea 
of the proof of Theorem 3.10 is to show that the sequence (p k )k£N is a Cauchy sequence in a 
complete metric space, hence converges. We use estimates between approximations p k and p m , 
to ^ k. Similar estimates are employed to obtain the result of Theorem 3.12. To prove Corollary 
3.11, we show that a solution at time t ^ 0 is provided by Theorem 3.10, if T > 0 is chosen such 
that t £ [0,T]. Moreover, this solution at time t is independent of the exact choice of T. 

4. Proofs of Theorems 3.10 and 3.12, and of Corollary 3.11. In this section we prove 
the main results of this paper: Theorem 3.10, Corollary 3.11 and Theorem 3.12. The essential 
part of the proof of Theorem 3.10 is provided by the following lemma: 


Lemma 4.1. For fixed po £ .A4 + ([0,l]) and (ak)ke N satisfying Assumption 3.5, the cor¬ 
responding sequence {p k )k^n defined by (3.3) is a Cauchy sequence in C([0,T];_A4 + ([0,1])). In 
particular, there is a constant C such that 


sup \\p k — pf ||bl ^ C 

tG[0 ,T] 


max 

jG{0,...,Y fc -l} 


(*i+i 



for all k, to £ N satisfying m ^ k. 


Proof. Fix k,m £ N with m ^ k, let t £ [0,T] be arbitrary and let j £ {0,... ,Nk — 1} be 
such that t £ (t k ,t k +1 \. Define, for appropriate ^ 1 ; the ordered set 

{n ■■ N&} := {t k } U (a m n (t k ,t k +1 ]^j U {t k +l }. (4.1) 

The set a m D (t k , t k +l \ contains all t™, £ £ {1,..., N m }, such that t k < t™ ^ t k +1 . For the sake of 
being complete, we emphasize that any duplicate elements that might occur on the right-hand side 
of (4.1) are not ‘visible’ in the set on the left-hand side. Assume that i £ {0,..., N& — 1} is such 
that t £ (rj, Tj + i]. To simplify notation, we write vf := v[pf ] for all k £ N and £ £ {0,..., N^}. 
Define io £ {0,..., N m } to be the smallest index such that t™ ^ t k . 
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Case 1: t™ = t). In this case, there is a q G {0,..., N m — 1} such that r, : = Hence, 
k = ^°E k, and k* = QrU 4 < • 

We estimate 

Ilk -klk < WQtrM* -OIIbl + II (kE - K-kklllk 


< K - Kllk exp (J^Il (r - n) + ||/||bl (r - n) e'^ (*-«) 

IIKIItv exp (l (t — Ti) + II/IIbl (r-Ti)e i(r ~ Ti) ^- 


+ K'n - v. 


0 u i IIOO 


(r - Ti ) + ( t - r l ) 2 ||/|| oc e^ll- {T ~ Ti) ], (4.2) 


using Corollary 2.9 and Lemma 2.10. Here, L denotes rnin(|k |l , |i>/|l). In view of Lemma 3.4, 
we define R := ||/zo||tv • expfll/Hoo T). From Lemma 2.7 (with s = 0), and Parts (ii) and (iv) of 
Assumption 3.1 it follows that 


l 

l|t>0-«riloo < Mr (lK -M™IIbl + E IK -^-JIbl) 

i =1 

i 

< Mr IK - K Ilk + Mr E IIQrflE, - ke-i Ilk 

£=1 

i 

< Mr UK - Ilk + Mr E ^ (ll/lloo + iffl) e"'"~ T k - r,_r) 

f=i 

^ M r IlK - Ilk + M* i? (ll/lloo + K r ) e"'"- T k - To). (4.3) 

We combine (4.2) and (4.3), and use Part (iii) of Assumption 3.1 and the basic estimates r — Ti ^ 
Tj+i — Ti and t,;_|_ i — Ti (in suitable places) to obtain that 

Ilk - klk < exp (A, k +1 - Ti)) IK - Kllk 
+ ^2 k+i - n) ||k 0 — Kllk 

+ A 3 (r i+ i - t ; )(t ; - r 0 ) ( 4 . 4 ) 


for some positive constants Ai, A 2 and A 3 that depend on /, T and R, but not on i or j. This 
upper bound holds for all t e (T*,Ti+i]. 

Case ,2: t k < t™ and i = 0. Note that j / 0 and io / 0 must hold. We recall the notation 
vf := v\n *J for all k G N and £ £ {0,..., N^}. In this case, 


l4 = Qr°- 


u k 

T 0 R'To ) 


and 


where v := v[n7L ]. Similar to (4.2), we have 
*0-1 


k 1 = Qk 


TO 




m 

'To’ 


Ilk - klk < IIQr-roCko - a 


IBL 


+ l|(Qr°-T 


Qkjk 


* 

BL 


< like - K Ilk exp (jklL (T - To) + II/IIbl (t - r 0 ) 

+ ||k - vU |K||tv exp (l (t - T 0 ) + ||/||bl (t - r 0 ) e L ^~^ ■ 

■ (T - To) + (t - T 0 ) 2 ||/|| oo e ll/l|o ° (T_ro) 


(4.5) 
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where L = min(|uo| L , Ml)- We define R := MoIItv • exp(11/11oo T)\ cf. Lemma 3.4. The analogon 
of (4.3) is 


H ~ elloo < Mr (|K - - MMbl) 

= Mr IK - KIIbl + Mr II Q%_r ^ - ^llk 
< Mr ||/4„ - /C 0 IIbl + Mr R (||/1|oo + Kr) ell'H” T (r 0 - f), (4.6) 

with f := t™_ i- Together (4.5) and (4.6) yield 

11/4 - /CIIbl < [exp (Ai (rr - to)) + A 2 (n - To)] ||< - KHk 

+ A 3 (n - t 0 )(t 0 - f) (4.7) 


for the same positive constants A\, A 2 and A 3 as in (4.4). Here, we used Part (iii) of Assumption 
3.1 and the estimates r — To < n — tq and n — r 0 < T. The upper bound (4.7) holds for all 
tG(t 0 ,ti]. 

Case 3: tj < t™ and i ^ 1. In this case, tj < n < tj +1 and hence there is a q € {1,..., N m — 1} 
such that Ti = t We have 


/4 = Qt-t- / 4 *> 


and 




Estimate (4.2) also holds in this case. Because t™ > tj there is no q G {0,..., N m — 1} such that 
tq = t™, and therefore u[-] is not to be evaluated at /x™. Consequently, we have instead of (4.3), 


vz-v? 00 < 00+ «r*-M 


£iK l -«£- 1 iik 


£=2 


< 


1*4 - Mloo + Mr || Q\_ f /X™ M^IIbL + Mr £ II Qrf-r^ ~ M™-! Ilk 


with u := v[/xjm J and r := t™_!- Note that the sum on the right-hand side might be empty. 
Using the idea of (4.3) and the result of (4.6), we obtain 

11*4 - i>nu < Mr UK - KIIbl + Mr R (||/||oo + K r ) e^ T (r 0 - f) 

+ Mr R (ll/IU + K r ) e'l'll” T (n - f) 

< Mr \\i4o - /C„ Ilk + 2 MrR (ll/lloo + A'M ell'll” T (r* - f )• (4-8) 

Due to (4.2) and (4.8), we have 

11/4 - /CIIbl < exp (A, (r <+1 - r,)) ||/x* - <||k 

+ ^2 (u+i — n) ||/4 0 - K Ilk 

+ 2 A 3 (Tj+1 - Ti)(ri - f) (4.9) 


for all r G (TjjTi+i], where Ai, A 2 and A 3 are the same constants as in (4.4) and (4.7). 

We now combine the estimates obtained in Cases 1, 2 and 3: it follows from (4.4), (4.7) and 
(4.9) that 
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for all i G {1,..., — 1}, while for i = 0 


sup ||/4 - /CIIbl < [exp {Ax (ri - r 0 )) + A 2 (n - r 0 )] ||/4 0 - /4 Jbl 

re(r 0 ,ri] 

+ A 3 #) (n-r 0 ). 


We have used that n - t 0 ^ M (k) in (4.4), r 0 - f < t% - t£_ x < M( m ) < in (4.7) and 
r, - f < t nU) - t 0 + t 0 - f ^ +t%- C-i < M (fe) + M (m) < 2 

in (4.9). This is the place where we use explicitly that partition a m is ‘finer’ (or: ‘not coarser’) 
than ak in the sense that M ^ ; cf. Assumption 3.5. By an induction argument one can 

show that the upper bound 


sup ||/4 - /CIIbl < II exp ( Al ( r « +1 “ T ?)) 

re(T 4 > r 4+1 ] U = ^ +1 


A 2 (r m - r/) H/4, - /£ ||k + 4 A 3 M< fc > (r /+1 - r<) 


II CXp ( Al ( T «+ 1 “ T <?)) Il^ro - llBL 
KQ—O ) 


(4.10) 


holds for all i € {0,..., — 1}. The products in brackets are equal to exp {A\ (r i+ 1 —Te + 1 )) and 

exp {A\ (ri+i — To)), respectively. By using these explicit expressions and by taking the supremum 
over i on the left-hand and right-hand sides of (4.10), we obtain 


iv°>_ l 


sup ||m4 — M^IIbl < Y ex P {Ax {t n u) — Ee+i))- 

' r 6(TO,T Jv0) ] 


i=0 


• [A 2 (r, +1 - r<) H/4, - /X™ ||| L + 4 A 3 M<*> (r m - r,)] 

+ exp (Ai (rjvO) - T 0 )) ||/4 0 - n™ 0 ||bl- ( 4 - n ) 

Since t n w — r^+i < r N uy - t 0 for all £ £ {0,..., — 1}, and r 0 = t k and t n u) = 4+i> it follows 

from (4.11) that 


sup Wi4 - (J?\\bi < 

T€(t k ,t k ,,] 
v 3 ’ J + 1 J 


N (j) -1 


1 + A 2 £ {n+x Tt) ■ exp (A, (4 +1 - 4)) ■ ||/4* - mJSIIbl 




jv<n_i 


+ 4 A 3 Af (fc) (4+i - 4) exp ( Al (^+i “ 4)) H ( r ^+i - r 4 

e=o 

= [1 + A 2 (4 +1 - 4)] • exp (Ai (4 +1 - 4)) • ||/4* - /i^llk 
+ 4 A 3 M <*> (4 +1 - 4) exp {Ax (4 +1 - 4)). (4.12) 


Hence, we have that 


sup ||/4 - /CIIbl < [1 + ^2 (4 +1 - 4)] • exp {Ax (4 +1 - 4)) • sup ||/4 - /CIIbl 

Te(iJ,tJ +1 ] 


+ 4 A 3 Af« (4 +1 - 4) exp (Ai (4 +1 - 4)) 

for all j £ {1,..., Nk — 1}, and for j = 0 we have 

sup H/4 - /CIIbl < 4 A 3 M<*> (4 - 4) exp (A x (4 - t k 0 )), 

re(t k ,t k ] 


(4.13) 


(4.14) 
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because /A = p 0 = p™. By an induction argument similar to the one leading to (4.10), we obtain 
that 


j j j 

sup ||/4 -/CIIbl<E n [l + A 2 (t k q+ 1 -t k q )\-e^(A 1 (t k q+ 1 -t k q )) 

re (*M+l] fcO \9=^+l 

■ 4 A 3 M(‘) (f , fe +1 - i*) exp (A, (t k +1 - t k )) (4.15) 

for all j G {0 ,..Nk — 1}. Note that [l + A 2 (t k +1 — t k )] < exp ( A 2 (t k +1 — t k )) for all q G 

{0,..., Nk — 1}. Define A 4 := A\ + A 2 . It follows from (4.15) that 

l 

exp (A 4 (t k j+l - if +1 )) • (i£ +1 - tf) 

T)(t k + 1 -t k ). (4.16) 

We take the supremum over j on both sides of the inequality (4.16) and get 

sup ||/4 -/CIIb L < 4A 3 M«Texp((A 1 +A 4 )T). (4.17) 

re[0,T] 

Note that we extended the supremum from r G (0,T] to r G [0,T], but this does not change 

the upper bound. Define C := 4 A 3 T exp((A 3 + A 4 )T) to get the result of the lemma. Because 

7 Vf( fc ) —»• 0 as k —> 00 , the estimate (4.17) implies that (p k )ken is a Cauchy sequence. □ 

Remark 4.2. It is crucial that we use || ■ ||g L and not || ■ ||tv in Lemma 4-1. The factor 

11 Uq — v^Woo appears in (4.2) due to Lemma 2.10. Due to Assumption 3.1 (iv) we subsequently 

obtain an estimate in which \\p k 0 — H™ Hbl appears. Analogous estimates apply to ||uo — u||oo in 
(4.5). Note that Lemma 2.10 builds on Lemma 2.2. In (2.13) the Lipschitz property of the test 
functions is explicitly used and hence, there is no direct way to formulate the result of Lemma 2.2 
in terms of || - ||tv- Consequently, we do not have an estimate of \\p k — p k+1 \\TV against ||v —i/||oc 
comparable to (4.2). 

We are now ready to prove Theorem 3.10. 

Proof. By definition, A4([0,1]) is complete in the metric induced by the norm II ' IIbl- The 
space A4 + ([0,1]) is a closed subspace of A1([0,1]), so A4 + ([0,1]) is complete. Hence, the space 

{v G C([0,T]; A4 + ([0,1])) : !/(0) = po} 

is complete for the metric defined for all p, v G C([0, T\; A4 + ([0,1])) by (3.6). 

For each initial measure po G A4 + ([0,1]) and for each k G N, consider the Euler approxima¬ 
tion p k defined by (3.3) corresponding to partition a*. This approximation p k is an element of 
C([0, T\; A4 + ([0,1])), because the semigroup (Qt)t^ 0 preserves positivity for all v G BL([0,1]); see 
[EHM15a, Corollary 3.4]. In Lemma 4.1, we showed that for given (ak)ken the sequence (p k )ken 
is a Cauchy sequence in {v G C([0, T]; A4 + ([0,1])) : u( 0) = po}, which is a complete space, as was 
argued above. Hence, the sequence ( p k )k <=n converges in {v G C([0, T\; A4 + ([0,1])) : i/(0) = ^o}- 

The limit is independent of the sequence of partitions chosen from the class characterized by 
Assumption 3.5. If (ctk)k£N and (/3k)ke N are two such sequences, then it is possible to construct 
a sequence ( 7 fc)fceN that has a subsequence that is also a subsequence of (afc)fc 6 N, and that has 
(another) subsequence that is a subsequence of (j3k)k&i- Moreover, ( 7 fc)fc<=N can be constructed 
such that the corresponding sequence of maximal interval lengths is nondecreasing. 



sup II p k - p ™|| <4 MMM exp(A 3 T) £ 

’’«#$.*?+ il 


e=o 

3 


<4 A 3 M^ exp(AiT) ^ 


e=o 


<4 A 3 M (fe ) exp((Ai + A 4 ) 
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Let (/« Q ' t )fc£Ni G u ' 8fe )fce n and (/t 7fe )fceN be the corresponding sequences of Euler approximations. 
The sequence (/z 7fc )fce^ can be shown to converge to the same limit as (/z Qfe )fc g N, and to the same 
limit as ([i^ k )km- Hence, ([i ak )ken and (k'OfceN converge to the same limit. This finishes the 
proof. □ 

The proof of Corollary 3.11 builds on the result of Theorem 3.10. 

Proof. Fix t ^ 0 and let T > 0 be such that t G [0,T]. For given /io G Al + ([0,1]) and v : 
A4([0,1]) x [0,1] —> R satisfying Assumption 3.1, a unique mild solution [i G C([0, T]; A4 + ([0>1])) 
exists, hence fit, the solution at time t, exists. We now show that this [it is independent of the 
choice of T. 

Let Ti, T 2 > 0 and assume without loss of generality that Ti <T 2 . For the given [i 0 G Ad + ([0,1]) 
and v : A4([0,1]) x [0,1] -A R, let [i denote the mild solution in CQO,Ti]; Ad + ([0,1])) obtained 
by partitioning [0,Ti]. Take a sequence of partitions (afc)fceN C [0,Ti] satisfying Assumption 
3.5, with corresponding Euler approximations (fi k )kew- Next, construct a sequence of partitions 
(/?fc)fceN C [0,T 2 ] satisfying Assumption 3.5, such that ctk C (3k for each k G N. More specifically, 
restricted to [0, Ti] each partition /3k coincides with ak- Note that such ((3k )fceN exists. Let (n k )ken 
be the sequence of Euler approximations corresponding to (/3k) ke N- 

For each k G N, the restriction v k | T ^ is defined by (3.3) with respect to the partition /3fcfl[0, T\] = 

ak- Hence zx fc | ^ is defined in the same way as /i fc , and thus sup Tg [ 0 Ti j \\[i k — v k Hk = 0 or simply 

u k \ T , = [if. Consequently, the same must hold in the limit as k —> oo, because of the triangle 
inequality: 

sup II^T -WIIbL < SU P \\Pt- U t\\bL+ sup ||Mr-Mr|lBL+ SU P II^-^tHbl 

t€[0,Ti] t€[0,Ti] re[0,Ti] t£[0,Ti] 

< sup 11/4 - ^tIIbL + sup ||^-Mt||b L + sup II^-^tIIbL- 

re[0,Ti] tG[0,Ti] t£[0,T 2 ] 

'-V-' '---V-' '-V-" 

=0 ->0 -X) 

So, sup Tg [ 0 Tl ] \\[i T — ^t||bl = 0- Hence, [i T = v T for all r G [0,Ti] and thus the solution at time r 
is independent of the final time chosen. □ 

Finally, we prove Theorem 3.12: 

Proof Let the mild solutions /t and v be given and let (ccfc)fceN be an arbitrary sequence 
of partitions of [0, T] satisfying Assumption 3.5. Let (/ i k )km and (i^ k )keN denote the sequences 
of Euler approximations defined by (3.3), both for the sequence of partitions (ctk)kefh and with 
initial conditions [ 1 q and z/o, respectively. 

Since /z and v are mild solutions 

/i = lim [i , and 

k—y oo 

v = lim v k 

k—> oo 

hold, with convergence in the metric (3.6). It follows from Lemma 3.4 that all elements of 

{[i k : k G N, t G [0, T]} U {v k : k G N, t G [0, T]} 

are bounded by R := R expdl/Hoo T) in both || ■ ||tv and || • ||g L . Fix k G N, let ak := {/§,..., t k Nk } 
and take j such that r G (t k ,t k +1 \. Consider the difference \\[i k — v k ||g L - 

We use an estimate in the spirit of (4.2)-(4.3)-(4.4). Note that the proof of Lemma 4.1 also 
holds if k = m, which implies = 1 and hence i = 0. It follows from (4.2)-(4.3), with i = 0 
and with u k instead of [i m , that 

11/4 - *4|Ibl < [1 +B 2 (t k + 1 - t k )\ exp (b, (t k + 1 - t k )^j ||/4j - u k k Ilk 


(4.18) 
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for some positive constants B i and B 2 that depend on /, T and R , but not on j or k. This 
estimate holds for all r G and resembles (4.4). We take the supremum over r G (tf,£f +1 ] 

on the left-hand side of (4.18), apply this relation recursively and take the supremum over j to 
obtain that 


( JV fe —1 \ 

n [ 1 + b 2 ( t e+ 1 - $)] ex P (#1 (*f+i - if)) ) IIMo - ^oIIbl 

e=o J 

< ( II exp ( B 2 (*f+i _ if)) ex P ( B 1 (if+i - if)) j IImo - J'oIIbl 

< exp ((Hi + B 2 ) (t% k - if)) \\p 0 - ^oIIbl 

= exp((B 1 + B 2 ) T) \\ho - ^IIbl, (4.19) 

for all k G N. The triangle inequality yields 

sup ||/Z t -I/ t ||bl < SU P ll^r “ ^tIIbL + SU P IIMt-Mt||bL+ S U P ||^-^t||bl. 

r£[0,T] t6[0,T] t£[0,T] t£[0,T] 


->0 


whence the same estimate as in (4.19) holds for sup Tg [ 0 T , \\/j, T — p t |Ibl- ^ 


Remark 4.3. We would have been inclined to use directly (2.20) on the interval (tj,tj +1 ], 
instead of deriving (4.18). We need, however, the exact dependence on (tj +i —tj) of the prefactor, 
to make sure that - after iteration over j - the prefactor in (4.19) is bounded. This dependence is 
not (directly) provided by (2.20), nor by the proof of [EHM15a, Proposition 3.5]. 


Remark 4.4. The result of Theorem 3.12 relies - via Corollary 2.9 and Lemma 2.10 - on 
Gronwall’s Inequality. This is possible here because we restrict ourselves to Lipschitz perturbations. 
In our previous work [EHM15a] we considered the more general class of piecewise bounded Lipschitz 
perturbations. Hence, there we stated explicitly (see the paragraph before [EHM15a , Proposition 
3.5]) that the standard approach did not work. 

5. Discussion. In this paper we have generalized the results of [EHM15a] to measure- 
dependent velocity fields via a forward-Euler-like approach. Our motivation was to derive flux 
boundary conditions for situations in which the dynamics are driven by interactions. Such dy¬ 
namics are in general more interesting than the dynamics that follow from prescribed velocity fields 
as in [EHM15a]. We managed to obtain a converging procedure, but only for bounded Lipschitz 
continuous right-hand sides. Hence, compared to [EHM15a], our results hold e.g. for boundary 
layers in which mass decays, but not for the limit case of vanishing boundary layer. We start off 
this discussion section (see §5.1) by commenting on the possibility to extend to piecewise bounded 
Lipschitz right-hand sides and to obtain the limit of vanishing boundary layer. Secondly, we point 
out (in §5.2) how this paper generalizes the results of [Evel5, Chapter 5] and how a number of 
open problems mentioned in [Evel5, Section 5.5] are now resolved. Ultimately, we suggest possible 
future research (§5.3). 

5.1. Piecewise bounded Lipschitz perturbations. To obtain the technical results in 
§2.3, we explicitly used the assumption that the perturbation / is bounded Lipschitz on [0,1]. 
Theorem 3.10 and Theorem 3.12 rely on the results in §2.3. We would have liked to obtain these 
results for piecewise bounded Lipschitz /, in particular to model decay of mass at one of the 
boundaries only (cf. [EHM15a]). In [EHM15a] we circumvent the arising problems by providing 
the solution explicitly in [EHM15a, Proposition 3.3]. In the setting of the present paper, this 
explicit form would be given for each interval (£*,i* +1 ], k G N, in (3.3) by 
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where v k := v[/i k k ]. In [EHM15a] we showed that it is possible to obtain the estimates needed 

to establish continuous dependence on initial data, because this explicit form has a regularizing 
effect on / and its discontinuities due to the integration in time. The key ingredient there, which 
is absent in the approach of the present work, is the fact that the velocity field is the same for all 
time. If one wants to prove Theorem 3.12 using (5.1) instead of the properties of the semigroup 
Q, one encounters that at some point for any At > 0 fixed a Lipschitz estimate of the form 

At At 

II J /($:(•)) ds-J /($“(•)) dslloo c ||u - v||oo (5.2) 

0 0 

is required, for all u and v taken from a class of admissible velocity fields. One would then proceed 
to estimate || u — i/loo against the BL-distance of the corresponding measures, using Part (iv) of 
Assumption 3.1. 

In view of [EHM15a], the restriction that the velocity should not be zero at discontinuities of 
/ is reasonable, but even if we are willing to obey that condition, an estimate like (5.2) cannot be 
expected to hold. Let f(x) = 0 if x £ [0,1) and /(1) = —1. Take e > 0 and take v = e, u = —e. 
Then (for e < 1/At) 

At At At At 

II J /($;(•)) ds-J /(<&“(•)) ds IU > I J /(<&:( 1)) ds-J /(*“( 1)) da\ 

0 0 0 0 

At At 

= I J JX 1 ) ds- j f( 1 - £ s ) ds| = At. 

0 0 

Since At > 0 is fixed and ||u — i/loo = 2e can be made arbitrarily small, (5.2) cannot be satisfied. 
An additional difficulty is that it remains to be seen how we can assure that a condition like 
u(l) ^ 0 is satisfied by a velocity field that depends on the solution itself. 

5.2. Uniqueness of mild solutions and generality of partitions. In [Evel5, Section 5.5] 
we point out that there are two reasons why we obtained uniqueness of mild solutions there. On 
the one hand, this is because the constructed approximating sequence converges, thus inevitably 
each subsequence (cf. Definition 3.8) converges to the same limit. This statement still holds true 
for the present work. On the other hand, uniqueness holds in [Evel5, Chapter 5] because there we 
only constructed one approximating sequence, namely by partitioning the interval [0,7/ into 2 k 
subintervals. In this respect, the present paper is a considerable improvement. The class of admis¬ 
sible partitions (see Assumption 3.5) includes partitions into q k equal subintervals for arbitrary 
q G N + ; see Example 3.6. We conjectured in [Eve 15, Section 5.5] that the sequence of correspond¬ 
ing Euler approximations converges, and the results of this paper confirm that conjecture. The 
fact that, in this case, each interval (t k , t k +1 \ is split into q m ~ k subintervals (t™, is generically 
treated by introducing the number and using a recursion over index i € {0 ,..., — 1} to 

obtain (4.11). In [Evel5, Chapter 5], however, we performed explicit calculations, using that each 
(t k ,t k + 1] is split into two subintervals. 

In [Evel5, Section 5.5] anticipated that using a sequence of non-uniform partitions of [0, T], would 
imply the need for a condition regularizing the variation in subinterval lengths to make sure that 
all subintervals become small sufficiently fast as k —> oo. In the present work we show that it 
suffices to have for the maximum subinterval length M^ —>• 0 as k —> oo. 

The iterative argument in [Evel5, Chapter 5] requires that the partition for index k + 1 is a re¬ 
finement of the partition for index k (more particularly: a division of each subinterval into two). 
The complications expected to occur if subsequent partitions are not refinements are resolved in 
the current work, by introduction of the index io in the proof of Lemma 4.1 and allowing for the 
case t k ^ t£. 
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The final contribution of the present work to be mentioned here is that in Theorem 3.10 we have 
positively answered the question posed in [Evel5, Section 5.5] whether the mild solutions obtained 
as limits of distinct sequences of partitions are actually identical. 

5.3. Future directions. The extension of the results stated in §3 to functions / with dis¬ 
continuities would clear the way for an approximation procedure like the one treated in [EHM15a]. 
That is, to have / nonzero only on the boundary of the domain and to approximate it with a se¬ 
quence of bounded Lipschitz functions (/ n ) n e n C BL([0,1]). In [EHM15a] we showed convergence 
of the corresponding solutions as n — > oo (for v £ BL([0,1]) fixed). The challenge would be (i): to 
establish the well-posedness of the problem for discontinuous /, and (ii): to show that the Euler 
approximation limit and the boundary layer limit commute. 

Let us focus on the vanishing boundary layer like in [EHM15a]. Assume there are regions around 0 
and 1 in which mass decays, and that these regions shrink to zero width. That is, there is a sequence 
(/n)ne n C BL([0,1]) and there is an / satisfying f(x) = 0 if x £ (0,1) and e.g. /(0) = /(1) = —1, 
such that /„ —> f pointwise, and the Lebesgue measure of the set {x £ [0,1] : f n {x) 7 ^ f{x)} tends 
to zero as n —> 00 . If we assume that we can extend the results of this paper to piecewise bounded 
Lipschitz /, then well-posedness for the limit case is guaranteed. It remains to be proven however 
that the solution for finite boundary layer actually converges to the solution of the limit problem. 
This is the same question as asking whether the two limits that we take, actually commute. The 
first limit is in the forward-Euler-like approach to obtain a mild solution. We assigned an in¬ 
dex k to the elements in the approximating sequence and proved in Theorem 3.10 that the limit 
“limfc_ ! . 00 ” exists (for / G BL([0,1])). The second limit “lim^oo” is the one involving the se¬ 
quence (f n )n£ N C BL([0,1]). Proving the well-posedness for / piecewise bounded Lipschitz, is the 
same as proving that the limit “limfc_>. 0 o lim n _»oo” exists. Proving that the sequence of solutions 
corresponding to each /„ actually converges to some limit in C([0, T]; yV( + ([0,1])) is equivalent to 
proving that “linin-Kx, lim^oo” exists. To conclude that the two limits commute, an additional 
argument is needed. It requires a characterization of “linin^oo limfc_ >00 ” that can be compared 
to “lim / c _ HX3 linin^oo”. Both proving that “lim^^oo lim^oo” exists and characterizing the limit 
can be a difficult task, however, since our current results do not provide an explicit expression for 
“linifc^oo”. A possible way to characterize the limit “lim^oo” could be to show that the mild 
solution obtained in this paper is actually a weak solution, and use the weak formulation of ( 1 . 2 ) 
as a characterization. If the solutions obtained in this paper are weak solutions, this is also a 
further justification of the terminology ‘mild solutions’. 

An additional result to be derived concerns the stability with respect to parameters, in par¬ 
ticular with respect to / and the specific form of v. Stability statements are essential in view of 
parameter identification. It is important to know how measurement errors in the parameters affect 
the solution of our model. In fact, Lemma 2.10 already provides stability in v for the solution of 
[EHM15a], provided that / G BL([0,1]). 

Moreover, we would like to study the long-term dynamics of the solutions 1 1—7 fit for various initial 
conditions. 
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